Preprint typeset in JHEP style - HYPER VERSION 



KEK-TH-1324 



LHC signature of supersymmetric models 
with non-universal sfermion masses 

Os ' 
O 
O . 

, Sung-Gi Kim 

^ , Department of Physics, Tohoku University, Sendai 980-8578, Japan 

o 

CO 



^1" 

m 

CN 



Nobuhiro Maekawa, Keiko I. Nagao 

Department of Physics, Nagoya University, Nagoya 464-8602, Japan 



*5 1 Mihoko M. Nojiri 

KEK Theory Center, IPNS, KEK, The Graduate University for Advanced 
D \ Studies (Sokendai), 1-1 Oho, Tsukuba, 305-0801, Japan 

IPMU, Tokyo University, Kashiwa, Chiba, 211-8568, Japan 



Kazuki Sakurai 

> 



KEK Theory Center, IPNS, KEK, 1-1 Oho, Tsukuba, 305-0801, Japan 

Abstract: We study the LHC signature of the minimal supersymmetric standard 
model with non-universal sfermion masses. In the model, soft masses of gauginos and 
the 3rd generation of 10 of SU(5) are around the weak scale, while other sfermion soft 
mass is universal and around a few TeV. Such sfermion mass spectrum is motivated 
not only from flavor, CP and naturalness constraints but also from E 6 grand unified 
model with non-Abelian horizontal (flavor) symmetry. The characteristic signature 
of the model at the LHC is the dominance of the events with 4 b partons in the final 
state together with high rate of mildly boosted top quark arising from gluino decay. 
The prominent high jet also arises from squark decay. We show it is possible to find 
the characteristic signature in the early stage of the LHC. The discrimination of our 
scenario from some CMSSM model points with similar signature may be possible 
with large integrated luminosity. The result of sparticle mass measurement using 
exclusive channel with the help of hemisphere analysis, and inclusive measurement 
of gluino and squark masses using Mt2 and M™ 2 m in some representative model points 
are presented. 
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1. Introduction 

Supersymmetry (SUSY) is one of the most promising candidates for physics beyond 
the Standard Model (SM). The minimal supersymmetric standard model (MSSM) 
has some attractive features, for instance improvement of the gauge coupling unifi- 
cation, the radiative electroweak symmetry breaking, and providing a dark matter 
candidate as the lightest superparticle (LSP) |], @, 0. Moreover, one of the most 
attractive features of the MSSM is the stabilization of the weak scale in the case that 
the SUSY breaking parameters and the higgsino mass parameter are around 
the weak scale. The signatures of the supersymmetry may be discovered and some 
properties of the MSSM will be revealed at ATLAS and CMS experiments at the 
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CERN Large Hadron Collider (LHC). The LHC signatures of various SUSY models, 
such as minimal supergravity (mSUGRA) model [||, |], ||, [?]], gauge mediated SUSY 
breaking (GMSB) model P, f?[ [T0[], anomaly mediated SUSY breaking (AMSB) model 
mixed modulas anomaly mediation (MMAM) model [|16|, |HJ , have 



been studied. 

The sparticle mass measurement is an important physics target at the LHC. 
Various methods have been developed for sparticle mass determination from event 
kinematics [4-7,18-35]. Endpoint methods of the various leptonic exclusive chan- 
nels is known to be very successful [§, |5], |(| ^ [18], By combining the measured 
endpoints of the invariant mass distributions of the jets and leptons from relatively 
clean and long cascade decay channels involving neutralinos and sleptons (/), 
one can determine not only the masses of the squark and gluino, but also the masses 
of neutralinos and sleptons arising from their cascade decays. Recently, progress 
has been made in the use of M<r 2 distributions for the sparticle mass measurement 
jZ7|, f?g, 0, 0, ^, g 0. It has been pointed out that the endpoint of the 
Mt2 distributions as a function of the test LSP mass (x) may exhibit a kink, which 
indicates the initially produced sparticle mass and the true LSP mass simultaneously 
Moreover if we define the My 2 and the sub- system M^ 2 [j33, 34] inclu- 
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sively, we can roughly measure the squark and gluino masses even in the very early 
stage at the LHC [E2|, [ 



36 



The purpose of this paper is investigating the LHC signature of SUSY models 
with non-universal sfermion masses. Most of the SUSY models that have been stud- 
ied so far have universality of sfermion soft masses in the flavor space except for 
some literature ^71 Introduction of the sfermion non-universality may in- 

duce unacceptably large flavor changing neutral currents (FCNCs) and electric dipole 
moments (EDMs) f|0], 41, ^2], 43]. However, the sfermion non-universality may be 
partially introduced for the 3rd generation sfermions because constraints from the 
3rd generation FCNCs are not so severe |44, ^5], (See also Refs. 
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If soft masses of gluino and right and left-handed stops are around the weak scale 
while the other sfermion masses are around a few TeV, some of FCNC and EDM con- 
straints are relaxed with keeping weak scale stabilization. In grand unified theories 
(GUTs) such as SU(5) GUT, both the left and right-handed stops are involved in the 
3rd generation of 10-plet of SU(5). In this paper, we consider a minimal non-univeral 
models in which only the 3rd generation sfermions involved in 10 have a different 
soft mass (m^o) from the other universal soft sfermion mass (mo) at the cutoff scale. 
The non-universal sfermion mass scenario is motivated both on phenomenological 
and theoretical grounds. There are models that predict the non-universal struc- 
ture adopted in this paper along with realistic fermion masses and mixing matrices 

MM- 

This paper is organized as follows. In Section 2, we introduce a non-universal 
sfermion mass scenario and explain motivations. We identify the motivated region 
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of parameter space from some phenomenological constraints. In Section 3, we inves- 
tigate the LHC signature at some representative model points. The charactaristic 
signature of the model is the 4 b partons in the final state and mildly boosted top 
quark with high rate for gluino gluino production. The prominent high jet also 
arises for squark gluino co-production. We demonstrate that it is possible to find 
the characteristic signature of the model through particle level Monte Carlo Simula- 
tion with detector smearing. In Section 4, we search constrained MSSM (CMSSM) 
parameter space ||52|| , where all scalar fields have common SUSY breaking mass at 
the cutoff scale, and find a model point whose signature is similar to that of our sce- 
nario. We propose a key measurement to discriminate our scenario from the CMSSM 
model point. In Section 5, we discuss sparticle mass measurement in our scenario. 
We study leptonic exclusive analysis, the top reconstruction from gluino decay, and 
measurement of the tb endpoint. We also study inclusive Mj-2 and M™ 2 m distributions 
[ |3^ , |33| , |36| and demonstrate that gluino and the first two generation squark masses 
can be measured in a stage of the LHC such as J Cdt = 5 — 20 fb -1 . Section 6 is 
devoted to the conclusion. 



2. Modified universal sfermion mass scenario 

In supersymmetric models there is no quadratic divergence in the Higgs sector. 
Therefore the naturalness problem in the SM is significantly relaxed. However the 
supersymmetry is softly broken and the scale of the quantum correction to the Higgs 
mass is of the order of the soft SUSY breaking parameters. If the SUSY breaking 
scale is much lager than the weak scale, unnatural tuning among the soft masses and 
higgsino mass, /i, is required. 

Let us look this issue more closely. In the MSSM, the condition for the elec- 
troweak symmetry breaking (EWSB) is given as 



/r -<(A)-Am^ + 0(-^), (2.1) 



= -lul 2 



where m 2 Hu (A) is a soft mass of the up-type Higgs boson at the cutoff scale A, and 
Am|^ is a quantum correction to the rn 2 Hu at the weak scale, which is roughly given 

as 

. ? 61 Y t \ 2 9l /A 2 \ , . 

Ami i—f-m 2 - hi ( — k ) , 2.2 

Hu (An) 2 * Vmf/' V ; 

where Y t is the top Yukawa coupling and mj is the averaged stop mass. The Am 2 Hu 
is large if mj is much larger than the weak scale. In that relatively large 

cancellation is required among the terms in the right hand side of Eq. ( [2.1| ). 

Unlike stop masses, the other squark and slepton masses do not affect the Higgs 
potetial because their Yukawa couplings are small unless the bottom Yukawa coupling 
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is as large as the top Yukawa coupling. As long as both left and right-handed stop 
masses are around the weak scale, we can take their masses much larger than the 
weak scale. In SU(5) GUT, two stops are unified into a single 10-plet field at the 
GUT scale. We consider a model in which a soft mass of the 3rd generation of 10 
(IO3) is independent of the other universal sfermion soft masses at the GUT scale. 
We parameterize sfermion soft mass matrices at the cutoff scale as follows: 



■2 



m{ = ml , = ml . (2.3) 



5 

m 30/ V 



Here, 10 = (Q,U C ,E C ), 5 = (D C ,L). In this paper, we take the cutoff A to be the 
unification scale of the three gauge couplings (A~2 x 10 16 GeV). 

In the super-CKM basis, sfermion mass matrices for 10 sector are given as m~ = 

V}ml V f (f = u L ,d L , u R , e R ) at the cutoffs scale. Here V f are unitary matrices that 
diagonalize the Yukawa matrices as V^ L Y U V* R = Y^ mg . The mixing induced by Vf 
should be sufficiently small to avoid large FCNCs. Thus in addition to Eq. ( |2.3| ), we 
assume |7|, M M, 



{V UL )ij, (Vdjij, (Kjij, {V eR )ij < (VhKM)ij, {i^j). (2.4) 



In E 6 SUSY GUT models with SU(2) x U(l) horizontal symmetry H, Eqs. 
and (|2.4|) are derived along with realistic fermion masses and mixing matrices |E0| 



5jJ . In the models, all the 5-plets of three families and the first two generations 
of 10-plets (10i and 10 2 ) are involved in a SU(2) H doublet, ^(27, 2) a = (^1,^2), 
while 10 3 and MSSM Higgs fields are involved in SU(2) H singlets, ^(27, 1) 3 and 
$(27,1), respectively. 1 The Yukawa terms for \^(27, 2) a are forbidden by SUSY 
and the horizontal symmetry of the superpotential. On the other hand, a term 
W 3 ^(27, 1)3^(27, 1) 3 $(27, 1) is allowed. Thus, the Yukawa coupling for 10 3 can 
be of order one, which is identified as the top Yukawa coupling. 

Other Yukawa couplings arise though the higher dimensional operators after 
breaking the horizontal symmetry. They are suppressed by factor ((F) /A) n , where 
(F) is a breaking scale of H and n is some integer. Thus, the Yukawa hierarchy in 
10 sector is larger than that in 5 residing in ^(27, 2). This explains why the mass 
hierarchy of the up-quark sector (Q, U c C 10) is the largest and the lepton flavor 
mixing (L C 5) is larger than the quark flavor mixing. 

The model predicts the non-universal sfermion masses Eq. (|2.3|) at the leading 
order. The renormalizable sfermion soft mass terms can be written as 

= ^00(27, 2)t*0(27, 2) a + 77^10(27, 1) 3 | 2 , (2.5) 



lr The numbers in the parenthesis denote representations under E 6 x SU(2)h 
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where <fi represents the scalar components of the superfield Again, the model has 
not only a partial (lst-2nd) universality of 10-plet sfermions but also a full (1st- 
3rd) universality of 5-plet sfermions. The full universality of 5 sfermions is crucial to 
satisfy the FCNC and EDM constraints because the unitary matrices that diagonalize 
the Yukawa matrices for 5s are expected to have large off-diagonal entries like the 
MNS (Maki-Nakagawa Sakita) matrix. 

In the following, we identify a region of parameter space where both the natural- 
ness and the FCNC constraints are satisfied. In this model, the condition Eq. (|2.1| ) 
can be expressed numerically in terms of the fundamental parameters defined at the 
cutoff scale as [^3] 



774 ~ -1.9|/i(A)| 2 - 1.2m^(A) + 1.5m| + 5.9m 2 /2 + • ■ ■ . (2.6) 

The term proportional to min arises through the stop mass dependence of the RGE. 
The large cancellation is not required in the right hand side of Eq. ( [2.6|) as long as 
parameters /z, m Hu , m 30 and mi/ 2 are around the Z boson mass scale. Therefore we 
study the region of the parameter space where 

mi/a < 300 GeV, m 30 , m Hu (A) < 500 GeV. (2.7) 

In this case, \i is typically 200 to 500 GeV. 

The other sfermion masses are given by mo- If mo is also around the Z boson 
mass scale, FCNCs and EDMs severely constrain the flavor off-diagonal terms of 
the sfermion mass matrices and CP violating phases of the various SUSY breaking 
parameters. Such constraints are sometimes problematic to construct explicit models, 
because various sources that violate the universality are expected. 2 On the other 
hand, if m is much larger than the weak scale, the constraints can be relaxed. 3 

However, there is upper bound on itlq. A large mass splitting between the first 
two generations and the 3rd generation sfermions tends to make the 3rd generation 
sfermion mass squared negative through the 2-loop RG effects, and cause the color 
and charge breaking (CCB) ^5|. In addition, large tan/5 (tan/3 = (H u ) / (H d )) 
potentially cause the CCB problem in our scenario. The 10 3 couple to the 5 3 through 
the bottom Yukawa coupling, Y b . The negative correction to (m~) 33 is roughly given 
as 

AK)33~-g«il»(^> (2-8) 

Note = mo at the cutoff scale. This contribution would be as large as the 2-loop 
RG effect unless Y b is sufficiently small. Small tan /3 (tan f3 ~ 10) is also preferable in 



2 The RG effect or the effect of the gravity mediation is one of the sources to produce the 
non-universality. In models with the horizontal symmetry, the effect of the horizontal symmetry 
breaking is also the source to produce the non-universality. 

3 The constraint from the up-quark (C)EDM is still severe because some contributions do not 
decouple with increading mo . A spontaneous CP violation mechanism may solve this issue p9[ . 
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Figure 1: The allowed region on (m,x/2 — mo) plane. The other parameters are taken 77130 = 
m Hu (A) = m Hd (A) = 300 GeV, A = -600 GeV, tan/3 = 10, sgn(/j) = +. 



Eq SUSY GUT model with horizontal symmetry. In the models the bottom Yukawa 
coupling, which is originated from the Yukawa couplings for \&(27, 2), is suppressed 
because it is forbidden under the horizontal symmetry. 

In Fig. [5, we show the allowed region on (m^ — mo) plane. Here, the low energy 



particle spectra are calculated using ISA JET 7.75 |5q| . We fix the other parameters as 
m 30 = m Hu (A) = m Hd (A) = 300 GeV, A = -600 GeV, tan/? = 10 and sgn(/j) = +, 
where A is the universal trilinear coupling. In the black region the lighter mass 
eigenstate of stops, £1, is unacceptably light due to the 1 and 2-loop RG effects. 
We find that mo cannot exceed 2 TeV in the region where m.30, mi/2 ^ 300 GeV and 
— A > 600 GeV. Thus, in this paper we study the parameter region where 4 

m ~ 1 - 2 TeV, tan (3 = 10. (2.9) 



In the following, we call the scenario characterized by Eqs. (|2.3|) , (|2.7f) and (|2[ 
modified universal sfermion mass (MUSM) scenario. 

The mass of the lightest CP-even Higgs boson, m^, tends to be lighter than LEP 
II SM Higgs mass bound m VSM > 114.4 GeV [|57|]. To push up above the bound, 
the quantum correction to the Higgs quartic coupling is crucial p9 . This typically 



requires large stop masses or a trilinear coupling. Since large 77130 and mi/2 are not 
preferable in view of the naturalness, we search an allowed region in the direction 
of large \A§\. Fig.[| shows the allowed region on (m^ — A ) plane. The red region 
is excluded by the condition m h > 114 GeV. We find that the Higgs mass bound 



4 This scenario however cannot explain the anomaly of the muon g — 2 J54j, because [i and 
are heavy due to the large m . 




180 200 220 240 260 280 300 320 340 



m 1/2 [GeV] 

Figure 2: The allowed region on (m 1 / 2 — 
Aq) plane. The other parameters are taken 
as mo = 1 TeV, TO30 = rriH u (A) = m,H d (A) = 
300 GeV, tan/3 = 10 and sgn(^) = +. 




180 200 220 240 260 280 300 320 340 



m 1/2 [GeV] 

Figure 3: The allowed region on (m 1( / 2 — 
Aq) plane in the light Higgs scenario. The 
other parameters are taken as mo = 1 TeV, 
m 30 = 200 GeV, tan/3 = 10, \i = 250 GeV 
and niA = 105 GeV. 



requires \A \ > 400 — 600 GeV in the region where m 30 , ~ 200 — 300 GeV. The 
two black dashed lines represent /i = 300 GeV and 500 GeV at the stop mass scale. 
The (jl value is not sensitive to A compared with its dependence on m 1 / 2 in Fig. 0. 
This means the weak scale mn u is not so sensitive to Aq, and we can take large value 
of A Q without making naturalness worse. However a large \A \ may also cause the 
CCB problem, because it may lead one of the stop masses squared to be negative. In 
the black region of Fig. 0, the lighter stop becomes unacceptably light, < 95 GeV 
pjfl , due to the large \A \ value. 

To perform Monte Carlo simulation studies we choose a representative parameter 
set, Point A with m = 1000 GeV, m 30 = 300 GeV and m 1/2 = 270 GeV. The other 
parameters are listed in Table [1]. In order to see the m dependence of the collider 
signature, we also choose model points AH1 and AH2 with m = 1400 and 1700 GeV, 
respectively. These model points are shown in Figs[l], |2| and Table [l]. Point U shown 
in Table|l| is a CMSSM model point which will be investigated in Section f|. The 
masses of some sparticles are shown in Table ^. 

If Higgs masses can be non-universal, there is another allowed region where the 
heaviest CP-even Higgs boson, H, becomes the SM-like Higgs boson |60], [61| . This 
scenario is referred as a light Higgs scenario or an inverted hierarchy scenario. This 
can be realized if the CP-odd Higgs boson mass, rriA, satisfies itla ~ 100 GeV, which 
also prefers small Because of the small tjia and \i parameter, the thermal relic 
density of the xl tends to be the same or small compared with the observed dark 
matter density p^L In this scenario all Higgs bosons, h, H, A, H ± , have small 
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Parameters 


Values 


A 


AH1 


AH2 


B 


U 




1000 


1400 


1700 

-L 1 \J\J 


1000 


1000 

-L \J\J \J 


'"oil 


300 


300 


300 


200 


1000 


mi/2 


270 


270 


270 


200 


270 


4) 


-600 


-600 


-600 





-1600 


tan/3 


10 


10 


10 


10 


20 


m 2 Hd (A) 


(300) 2 


(300) 2 


(300) 2 


-(216.0) 2 


(1000) 2 




300 


300 


300 


196.7 


1000 



Table 1: The parameters of each model point. The unit of mass parameters is GeV . 





Masses (GeV) 


Particles 


A 


AH1 


AH2 


B 


U 


q 


1150 


1500 


1780 


1080 


1145 




321 


262 


187 


296 


281 


h 


540 


499 


456 


400 


856 


9 


697 


711 


721 


537 


706 


X? 


110 


111 


111 


77 


114 



Table 2: The masses of some sparticles at each model point. 



masses around the weak scale, therefore they would contribute to various rare B decay 
processes. Actually, constraints from Br(B+ — > t + u t ) |54] and Br(B s — > p5| 
exclude tan/3 > 15 region pBfl . And Br{b — > 57) requires sign(/i) = + and small 
masses of both ti and xt [§3 , so that the charged Higgs - top quark contribution can 
cancel with the chargino - stop contribution. The prediction of Br(b — > 57) however 
depends on the off-diagonal entries of V UR , V UL and Vd L in Eq. ( [2.4| ) HIT 



In Fig.|^, we show the allowed region on (mi/2 — Aq) plane. Here we fix uiq = 
ITeV, m 30 = 200 GeV, tan/3 = 10, /i = 250 GeV and m A = 105 GeV. The heavier 
Higgs boson in the MSSM is the SM like in this case, and LEP II bound should be 
applied to H. The red region of Fig.|3| is excluded by the condition m# > 114 GeV. 
The allowed region is widely extended to the small A$ region compared with the 
normal case. We choose Point B defined in Table [I] as a representative parameter 
point. 
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3. The characteristic signatures of MUSM 



3.1 The number of b jets: 

In MUSM scenario, m is much larger than m^. Then, the gluino 2-body decay 
mode into the first two generation squarks, g — > qq, is closed. On the other hand, in 
a wide parameter region t\ and &i are lighter than gluino due to the RG running and 
left-right mixing effects even if m 30 > m^. Then, g — > t^it and 616 modes entirely 
dominate the gluino decay. Since a gluino is a flavor singlet, the gluino decay chain 
contains at least 2 b jets (6-6 pairs). Therefore g-g and q-g(q) production events have 
4 6 partons. This characteristic feature can be observed at the LHC by counting the 
number of b tagged jets. 

To simulate the LHC signature at the model points selected in Section^], we 
calculate the low energy particle spectra and the sparticle decay branching ratios by 
ISAJET. The SUSY events are generated by the parton shower Monte Carlo HERWIG 
[3?], HH , and the detector resolutions are simulated by AcerDET |59f . We assume the 



collider center of mass energy is y/s = 14 TeV. Throughout this paper, we adopt the 
following SUSY cuts to reduce the SM background: 

• Nit > 4, AC > 1, 

• M cS = Eti \Pt\ +$t> 500 GeV, 



• Ifjf > max{200GeV,0.2M eff }, 

where p^) is the transverse momentum of i-ih jet (p^) > p^) for i < j) and N^ lQQ ^ 
is the number of jets with p T > 50(100) GeV and \rj\ < 3. 

The distribution of the number of b jets at each model point is shown in Fig. [|. 
Here we count the b tagged jets with prib) > 50 GeV and \rj(b)\ < 2.5, and require 
no isolated lepton in the event. The number of SUSY events corresponds to lib -1 of 
the integrated luminosity. We assume the b tagging efficiency is 60 %. To compare 
the distributions with those of the benchmark model points with universal sfermion 
masses, the distributions for the SPSla and SPS2 model points [[F(J are also shown 
in Fig.|j. The peak of the distributions is at zero at SPSla and SPS2, while zero b 
tagged events are suppressed at Point A, AH1 and B. 

The suppression of no-b jet SUSY events indicates the following mass relation: 



win > rn„ > 



Indeed in the distributions of mSUGRA benchmark points SPSla— SPS9, this feature 
is not seen (See Appendix|BD. Even at SPS2, where (rrig, nig, m^, m^) = (796, 1560, 
963, 1301) GeV and the gluino branching ratio into 2b + X is about 70%, the peak 
is at zero. 
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Point A 



300 r 
• 250 



Number of b— tagged jets 



Point B 




Number of b— tagged jets 



Point AH 1 



T 225 
^ 200 



Number of b— tagged jets 



SP5 I a 




Number of b— tagged jets 



Point AH2 



Number of b-tagged jets 



SPS2 




Number of b-tagged jets 



Figure 4: The distribution of the number of b tagged jets per lfb 1 at each model point. The b 
tagging efficiency is assumed to be 60%. 



An exceptional case is Point AH2. At this point, the distribution peaks at zero, 
although the mass relation Eq. (ft.lQ is satisfied. At this point the ti is as light as top 
quark, = 187 GeV, due to the large m value (m = 1700 GeV), and the 2-body 
decay modes t\ — > Xx b an d t\ Xit are closed. Then the flavor violating 2-body 
decay mode t\ — > Xi c dominates the t\ decay [[n]]. In such case, the distribution 
of the number of b jets are not helpful for the model discriminations. It has been 
shown in [72 that such a light stop may be detected at the LHC in 7 + 1$? or j + I$t 



channel. 

3.2 The highest px jet: 

In MUSM scenario, there is large mass splitting between mo and mi/2, m.30, but m 
is bounded above by the CCB constraint. The 1st and 2nd generation squarks, q, 
may be produced enough to be seen at the LHC. The number of SUSY events after 
the standard SUSY cuts at J Cdt = lfb -1 are listed in Table [| 

The SUSY production is dominated by g-g and q-g production processes at Points 
A, AH1 and B. On the other hand, the fraction of q-g{q) production is relatively small 
(10 %) at Point AH2, due to the very large mo value (m = 1700 GeV). The fractions 
of t-t and b-b productions are also small at all model points. These events hardly 
survive the standard SUSY cut due to the small masses of ti and b\. 

Once a heavy squark is produced, it decays mainly into g + q. The quark jet 
is expected to have relatively large transverse momentum since the mass difference 
between g and q is large. The order of the transverse momentum is around rriq/2. 
This quark jet tends to be the highest pt jet in the event. 
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Production 
Processes 


Product 


ion Ratios (%) 


A 


AH1 


AH2 


B 


U 


99 


32 


50 


55 


47 


36 


qg 


43 


22 


9 


38 


45 


qq 


7 


2 


1 


3 


9 


Mi 


9 


14 


19 


5 


5 


6161 


2 


5 


9 


3 





others 


7 


6 


7 


4 


5 


Total Events 


1484 


984 


1096 


3468 


1677 



Table 3: The number of SUSY events after standard SUSY cut at J Cdt — 1 fb Here we use 
the same character for a particle and the antiparticles. 



Point A Point AH 1 Point AH2 




(GeV) (GeV) (GeV) 

The highest pT (not b— tagged) The highest pT (not b— tagged) The highest pT (not b— tagged) 



Point B SPS I a SPS2 




(GeV) _ _ (GeV) (GeV) 

The highest pT (not b— tagged) The highest pT (not b— tagged) The highest pT (not b-tagged) 



Figure 5: The px distribution of the non-b tagged highest pt jet in the event at each model 
points. 



We show the pt distribution of the highest pt jet at each model point in Fig. [|[ 
Here we require that the rapidity of the highest pt jet is less than 1.5 (|?7(j lst )| < 1-5) 
and the jet is not h tagged. The number of the generated events corresponds to 
J Cdt = lfb -1 . In the figures the yellow histograms represent the contribution of 
g-g and t-t productions, while the shaded histograms represent that of g-q and q-q 
productions. The pt distribution of squark productions tends to be harder than 
that of gluino pair production. This is the contribution of the hard jet from the 
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Figure 6: The pr distribution of the b tagged highest pr jet in the event at each model points. 



heavy squark decay. Indeed, as m Q increases (mo = 1400, 1700 GeV at Points AH1 
and AH2), the peak of the shaded distribution moves to the high energy side. At 
the same time, the number of the events coming from squark production decreases 
significantly. 

The difference between the yellow and shaded distribution indicates 

m§- - > maxjmj - - m t(b) , - m^o} , (3.2) 

under the mass relation Eq. (|3.1|) . 

In order to see the contributions from gluino pair production separately, we can 
use the distribution of the b tagged jets. If the highest jet is a b jet, the jet is not 
originated from the first two generation squark decay. The contribution from the 
heavy squark decay can be removed by requiring that the highest pt jet is a b jet. In 
Fig.||, we show the pr distribution of the b tagged highest pr jet. The heavy squark 
contributions (the shaded histograms) are significantly suppressed in the figures. The 
difference between the distributions of non-6 tagged jets and b tagged jets suggests 
the large mass splitting between g and q. 

For comparison, we have done the same analysis for the benchmark model points 
SPSla— SPS9 (See Appendix^). There is no clear difference between the "6 tagged" 
and "non-6 tagged" distributions except at SPS2. At SPS2 the universal scalar mass, 
mo = 1450 GeV, is much larger than the universal gaugino mass, m 1( / 2 = 300 GeV. 
So the decay of the 1st generation squarks produce the hard jet, while the gluino 
3-body decays cannot produce hard b jets. 
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parameters are chosen as mo = m,H u (A) = mjj d (A) = 1000 GeV, min = 270 GeV, sgn(^i) = +. 



4. "Look alike" in CMSSM 

The characteristic signatures of MUSM scenario in the number of b tagged jets and 
the p!p distributions do not eliminate the possibility of the universal (CMSSM type) 
boundary condition of soft masses at the cutoff scale, though they are good indication 
of MUSM scenario. In this section we study a region of parameter space in CMSSM 
whose signature is similar to that of MUSM scenario, and seek for the outstanding 
observable that is useful to distinguish MUSM scenario from the CMSSM parameter 
region. 

If one takes m ^> mi/ 2 in CMSSM, the pj? distributions of b tagged and non-6 
tagged jets is similar to that of MUSM scenario as we have seen at SPS2. When 
one fixes itlq and mi/2, CMSSM still has the other 2 free parameters Aq and tan/?. 
We scan (tan/? — A ) parameter space with fixing mo = 1 TeV, mi/2 = 270 GeV 
and sign(/x) = + in Fig. [7]. In a gray region, g — > tit mode is open and it entirely 
dominates the gluino decay. Therefore there will be 4 b partons in the final state as 
in MUSM scenario. We choose such a model point U as a representative parameter 
point. The parameters are listed in Table [I] and the sparticle masses are listed in 
Table |^. 

We show the number of b tagged jets and the p^ distributions at Point U in 
Fig.|8|. The distributions are similar to those of MUSM scenario. This means the 
analysis presented in the previous section cannot discriminate MUSM scenario from 
the gray region in CMSSM, though the soft masses are universal at the cutoff scale. 
Therefore more information is required to distinguish these two scenarios. 
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Figure 8: The distributions of the number of b jets (Left) and the pr of the highest pr jet for 
non-b tagged(Center) and b tagged(Right) jets at point U. 



There are some differences in gluino decay branching ratios. In CMSSM case, 
the gluino decay can be dominated only by g — > tit mode. If one takes a large tan (5 
along with a large \Aq\, g — > bib mode is open as shown in Fig. 0. However ti is 
much lighter than bi in this region, and Br(g — > b\b) cannot be significant. On the 
other hand, Br(g — ► bib) can be as large as Br(g — ► tit) in MUSM scenario, unless 
ti is extremely lighter than bi due to the large \Aq\ effect. Therefore we can regard 
a sizable Br(g — > bib) as an indication of MUSM scenario. 

Distinction of g — > bib decay from g — > tit decay is not so easy at these model 
points. A difficulty comes from similarity of final states of these two decay modes. 
The main decay chains of t% and bi are as follows. 

■ • (a) bi -> tiW^ ■■■(d) 

(b) -xff-X?W«t •••(e) 

■■(c) ->3S&-xSZ«6 ■■■(/) (4.1) 

Here we ignore bi — > xV 3 mode because the branching ratio is tiny due to the hy- 
percharge and U(l)y gauge coupling suppression if bi ~ bi and the Xi is bino-like. 
In MUSM scenario or in the gray region of CMSSM, xt -> Hvl) and x° 2 -> ll(vv) 
modes are not open due to a large rriQ. Therefore we assume Xi ~* )$W an d 
X2 ~^ X\Z modes are open and dominate Xi an d X2 decay. Except for the decay 
modes (c) and (/), p — > tt and g — > bib have the same final state 2W + 2b + x^. 

The decay modes (c) and (/) are useful although their branching ratios are less 
than 10% in a wide parameter region. The final states of the gluino decay via the 
decay modes (c) and (/) are 2W + 2b + Z + x° and 2b + Z + Xi, respectively. If one 
of the W bosons decays into jets, the number of associate jets is differ by 2 between 
these modes. 

We show the distributions of the number of the jets with > 50 GeV and 
\r)\ < 2.5, (A^j ets ) in events with a Z — > iflj candidate at Point A and U in Fig^. We 
require there are opposite sign same flavor (OSSF) lepton pairs with > 15 GeV 
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Figure 9: The distributions of the number of jets in the event with Z — > ifl- at Point A and U. 



and \r)\ < 3 whose invariant mass satisfy \mu— mz\ < 5 GeV. 5 To reduce an additional 
Z boson source from q — > %%q — > ^Zq, we adopt g-g selection cuts: 

• The highest pt jet is b tagged, 

• pfi < 300 GeV. 

The shaded histograms in these figures represent the contribution from the q-g and 
q-q production events. The contamination from the squark decay q — > X2Q ~ * Xi^q 
is negligible after the 'g-g selection cuts. 

At Point A the peak is at N^ ets = 4, while the peak is at N^ cts = 6 at Point U. 
This suggests that the Z boson comes from b\ decay at Point A, while it comes from 
ti decay at Point U. By this analysis, we can discriminate MUSM scenario from the 
CMSSM parameter region, though a large integrated luminosity may be required. 6 



5. The sparticle mass measurement 

5.1 Exclusive analyses 

5.1.1 Conditions to have OSSF lepton pair signature in MUSM scanario 

SUSY particle masses may be determined by measuring the kinematical endpoints 
of sparticle decay products. Especially, the distribution of the events with OSSF 
leptons coming from x°2 -> if If -> Xilflf are useful [|, |, § % |l|, ||. However, 



there are several reasons for this channel may not be available in MUSM scenario. 



5 If X2 ^ x\Z mode is not open, the distribution of two lepton invariant mass would show 
kinematical edge induced by X2 XxH^ mode as shown in Section 5.1.2. In this case, we should 

cd 

hi 



require m?, dsc — mu < 10 GeV instead of \mu — mz\ < 5 GeV. 



6 Note in the region above a green line in Fig.Q, the decay mode (c) in Eq. (4.1) is closed. In this 



region, Z boson does not come from g decay in CMSSM, so the discrimination from this region is 
easy. 
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In this scenario gluino entirely decays into the 3rd generation squarks, t\ and 
b\. From a phase space consideration Br(t\ — > x^o) 3> 2Br(t\ — > x^t) an d Br{t\ — > 
X\t) ^> Br(t\ — > x^t), where the factor 2 in the first relation comes from the Dirac 
nature of the chargino. Thus, t% — > xfb and ti — > x?£ dominate the stop decay. If 
the gluino decay is dominated by g — > tit mode, x° does not appear with sufficiently 
high rate in the gluino cascade decays. 

In this scenario the b\ branching ratio into x° ma y a ls° be small, though gluino 
can decay into b±b. Assuming = 0, X2 = W3 an d &i = ^l, the tree level formulae 
of the b\ — > y§b and b\ — > t\W decay widths are given by [fft| 



. ml 
32n V to? 

b 



r^^-^-^h--^ , (5.1) 



9*^-3/2, .( \ f m t 



tiW)*^?\' (x,,x w ) Ij-. * , , - -f , (5.2) 



327T \(to~ — to- ) 2 + afmf / \m 



.ill!,/ \ 1 1 brir 



where = to- /mf , xyy = m^/mi and \(x, y) = 1 + x 2 + y 2 — 2(x + y + xy) and 

at = {At — //cot/5). In Eq. (|5.2|) T(6i — ► tiW) has a enhancement factor (jn 2 /m 2 v ). 
Thus, 61 decays dominantly into t\W unless it is kinematically suppressed. 

In MUSM scenario, the first two generation sfermions are much heavier than the 
second lightest neutralino. Therefore the 2-body decay mode of x° m t° the first two 
generation sleptons x° ~^ ^t^f ls kinematically forbidden. The OSSF leptons may 
arise from the 3-body decay mode x° Xi^t^T from the off-shell Z boson exchange. 
If the 2-body decay mode X2 ~^ Xi^ is open, it dominates the X2 decay. 

In summary, the sparticle mass measurement using the exclusive analysis with 
OSSF leptons may work if following conditions are satisfied: 

(i) g — * t\t mode does not dominate the gluino decay. 

(ii) b\ — > t{W mode is kinematically suppressed. (5.3) 
(iii) X2 ~~ * XiZ mode is kinematically forbidden. 

The three conditions can be satisfied when the mass differences m-r — m, and m^o — 
m^o are small, and it can be realized for small myz and \A$\. If all conditions are 
satisfied, a gluino cascade chain, 

g ^ ft l6 (i) ^ x°2b {1) b {2) - x\b {l) b^lfl-, (5.4) 



can have enough branching fraction. Here subscripts of the b quarks in Eq. are 
used to distinguish two b quarks appear subsequently. 

The low energy mass spectra and the decay branching ratios of various sparticles 
at Points A and B are shown in Table [| in Appendix^. At Point A, the three 
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Figure 10: The invariant mass distributions of ifh . The background distributions estimated 
from OSOF leptons are subtracted. The number of events corresponds to / Cdt =10fb _1 . 



conditions in Eq. ( |5.3| ) are not satisfied. On the other hand, at point B all of the 
three conditions are satisfied due to the very low SUSY breaking scales, m 1( / 2 = 
m 30 = 200 GeV and A = GeV. 



5.1.2 The exclusive analysis using xt ~^ Xi^th channel at Point B 

At Point B, SUSY events contain the cascade decay chain (|5.4j ). The OSSF lep- 
ton pair from the decay chain is relatively clean signal j|, |5|, ||, 0, [18], pl|. The 
background distributions of fake OSSF leptons can be estimated by distributions of 
OSOF leptons. The kinematical maximum of the two lepton invariant mass coming 
from X2 ~ * xV-tK a ^ Point B is given by 



m, 



X-2 



Xi 



61 GeV. 



(5.5) 



The distribution is shown in Fig|10|. Throughout this section, leptons are required 
to satisfy pr > 15 GeV and \t)\ < 3. 

For the decay chain (|5.4|), the invariant mass distributions m^, Tn l + rb(1) , m l + l - b ( 2 ) 
and rn m + r may determine sparticle masses. Their kinematical maxima at Point B 

i i 

are given as 



rr? max 

m bb 
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336 GeV, 



(5.6) 
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296 GeV, 



m^* =m~ Q - m 7 o = 460 GeV. 
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(5.7) 
(5.8) 

(5.9) 
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Since there are typically 4 b partons in the SUSY events, the invariant mass 
distributions suffer from the combinatorial background. To reduce the background, 



we adopt hemisphere method [[76], [77|| . The method divides jets and leptons from 
cascade decay chains into two groups called hemispheres whose entries are much likely 
to originate from the same mother particle. The groups are defined by hemisphere 
momenta, 

j£L = I>! , \ p!1, = £p. W , (5-10) 

i i 

where pf~' and pf^ are jet or lepton momenta that satisfy 

d(jP<£L,V? ] ) < dijp^J?), d(p® mi ,pf) < dip^pV). (5.11) 
Here the function d is defined by 

d(tfL,Pk) = (E^L - lp£Ll oose tt ) - , (5-12) 



(i) 

where 6^ is an angle between p hemi and pk- We require that jets involved in hemi- 
sphere satisfy pr > 50 GeV and \r]\ < 2.5 to reduce the contamination of soft jets. To 
find the two groups of jets and leptons that satisfy Eq. ( |5.10|) , we first choose the high- 
est pt jet and the jet with the largest prAR as pi(seed) and p 2 (seed), where AR = 
\J A<p 2 + At] 2 is the angle difference between the jet and the highest jet. Next we 
group jets into two groups under the condition d{p^> (seed), p^p) < d(p^ (seed), pf') 
etc.. Then axis momenta are defined as pW(ax) = ^2pf\ p^(ax) = YIpP ■ For the 
axis momenta, new groups are defined so that d(p^\a,x) , pf~ ) < d(p^ 2 \a.x), pf^) etc.. 
The procedure is iterated several times so that assignment converges. 

Although the probability that hemisphere correctly reconstructs original cascade 
chains is not so high, this method has an advantage for endpoint analyses. In the 
algorithm, two objects whose momentum directions are roughly the same each other 
tend to be in the same hemisphere due to Eqs. (|5.11|) and (|5.12j) . Because of this 



property, any invariant mass distribution of jets and leptons in the same hemisphere 
tends to be lower than that in the different hemispheres. Therefore wrong com- 
binations whose invariant mass exceed the signal endpoint are removed with high 
probability if we take jet pairs in the same hemisphere. 



We show the invariant mass distributions of bb, Ijl^b and l^l~bb in Fig. [11]. Here 
all b jets and leptons are required to be in the same hemisphere. If there are more than 
one or two b jet candidates, we take high pt b jets. We also require m l + l7 < m^L 
to reduce background OSSF leptons from Z boson decay. The distributions have 
endpoints near the theoretical expected values shown in Eq. (|5.6| ) to ( [5.9| ). However 
bb distribution does not show expected sharp edge structure due to the hemisphere 
selection. 
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Figure 11: The invariant mass distributions of bb (Left), blfl~ (Center) and lfl~bb (Right). 
The background distributions estimated from OSOF leptons are subtracted. The corresponding 
luminosities are J Cdt =50, 40 and 50 fb -1 for the left, center and right distributions. 



At this stage, we have 4 measured values m™^, m^ ax , rn^L and m?.™ for 
4 unknown sparticle masses m^o, m^o, and rrig. By solving 4 equations, ( |5.5| ), 
, ( |5.8| ) and (|5.9|), we can get all sparticle masses appeared in the decay chain 
in the stage around 50 fb -1 at Point B. 
In addition, we can check our results by selecting the events near ma endpoint. 
The two lepton system from the 3-body decay x° ~~ > Xi^C mus t be at rest in the Xs> 
rest frame when va l + l - = m°^i . The x% momentum can be estimated from a velocity 
of the two lepton system for the events with m l + l - < rn^L, if m^o is known ||, |T9[ |. 
Using the observed x° m ass from 4 endpoint measurements, (|5.5|), (|5lf), (|5.8|) and 
( |5.9| ), we can calculate the invariant mass distribution of x° + b from estimated x° 
momentum. For the decay chain (|5.4| ), the peak of the m^o b (2) distribution gives the 
b\ mass. We show the fn^o b distribution in Fig. [12] (Left). Here we use events that 
satisfy GeV < m™^t — m l + r < 10 GeV. Here we use all b tagged jets and leptons 

i i i i 

as the candidates because only a few events survive under the requirement that all 
b jets and leptons are in the same hemisphere. The distribution has a peak near the 
correct &i mass, m bi = 400 GeV. 

By using events around the b\ peak, we can subsequently estimate the b\ momen- 



tum, and calculate the invariant mass of b±b. For the decay chain ( |5.4j ), the g mass can 
be measured from the peak of the m b ^ distribution. We show the m bib distribution 
in the right figure in Fig. [12]. Here we use events that satisfy \m^ b k — m^o b \ < 15 GeV. 
The distribution has a peak near the correct g mass, m-g = 537 GeV. 

5.1.3 Top reconstruction and the tb endpoint at Point A 

The sparticle mass measurement is challenging for the class of the points A, AH1, 
AH2. The gluino entirely decays into t\ and b±, and decay modes of the t\ and b\ 



are dominated by (a), (b) and (d), (e) in Eq. ( |4.1| ), respectively. All these cascade 
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Figure 12: The invariant mass distributions of X2& (Left) and b\b (Right). In the left figure 
4 momenta of the \2 are estimated from the velocity of the lepton pair system whose invariant 
mass satisfies OGeV < — m ; + ; - < lOGeV by assuming correct X2 mass. In the right figure, 

momentum of y^b system whose invariant mass satisfies |to~o 6 — m x%b\ ^ 15 GeV is regarded as b\ 
momentum. The corresponding luminosities are J Cdt =30 fb -1 (Left) and 50 fb -1 (Right). 



chains have the same decay products 2b + 2W + xl, and a gluino pair system leads 
to 46 + 4W + 1$t- 7 Therefore, any exclusive analysis using jets suffers from a large 
combinatorial background. The gluino decay products contain a top quark with high 
probability. If the top quark is detected with a significant rate in SUSY events, it 
indicates the existence of the light 3rd generation squarks. 

The lighter stop mass is relatively light at Point A, mj =321 GeV, and the 
phase space of the decay is not small, rrig — (mr + m t ) ~ 200 GeV. The top quark 
from the gluino decay, g — > tit, is boosted, and the top decay products, b + 2j, tend 
to go in the same direction. The jets that goes in the same direction are efficiently 
picked up by the hemisphere method. In the analysis therefore we require all b + 2j 
are in the same hemisphere. 

Since there are a large number of jets in SUSY events, the combinatorial back- 
ground for W reconstruction still remains after the restriction. In order to estimate 
the background distribution from fake jet pairs whose invariant mass is ~ mw, we 
apply sideband subtractions [[[£], |7^] . We first define the W mass region and sideband 



regions as follows, 

\rrijj — m w \ < 10 GeV • • • (W mass region), 
\rrijj - (m w - 20 GeV)| < 10 GeV • • • (W sideband region I), 
\rrijj - (m w + 20 GeV)| < 10 GeV • • • (W sideband region II). (5.13) 

If the event contains several jet pairs in the same region, we choose the jet pair 



7 This event topology is essentially the same as in the focus point like region [^|. However in 
our scenario the branching ratio to the 46 + AW + I$t final states is very large because gluino can 
decay to the 3rd generation squarks in 2-body way. 
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Figure 13: Left; The invariant mass distribution of the jet pairs that give the closest value to 
rrijj = mw — 20 GeV for sideband region I, rrijj — mw for W mass region, rrijj = mw + 20 GeV for 
sideband region II. Center; The invariant mass distributions of jjb, where the jet pair is in the W 
mass region. The shaded region represents the estimated background distribution from the fake W 
jet pairs. The background is estimated from events in the sideband region I and II by rescaling the 
momenta and normalizing the number of events events. Right; The jjb invariant mass distribution 
after the background subtraction. 



whose invariant mass is the closest to mw for W mass region and to mw T 20 GeV 
for sideband region I (II), respectively. 8 In Fig. |TB| (left) we show the selected two 
jets invariant mass distributions. The number of events in W mass region is clearly 
bigger than those in the sideband regions. 

The central figure shows the jjb invariant mass distributions, where the two jets 
are the closest jet pair to mw- Here the open histogram shows the distribution of 
the events in W mass region. On the other hand, the shaded distribution shows the 
background estimated from the sideband events. Namely, the jet pair momentum in 
region I (II) is rescaled by a factor (mw/{m>w T 20 GeV)) before calculating rrijjb, 
and the two distributions are averaged. The right figure shows the W + b invariant 
mass distribution after subtracting the sideband distribution. It has a clear peak 
at the top mass, which indicates that the SUSY events contain top quarks with a 
significant rate. 

Because we find very prominent top quarks, it is natural to think tb distribution 
also show the clear kinematical structure. At pont A, the dominant gluino cascade 
decay chain is g — ► tit — > xfbt with 47% of the branching ratio. The m tb distribution 
should have an edge. The kinematical maximum of the tb invariant mass is given as 



9 9 \ 1/2 

m~ — mz± x 



mS 8 * = I mt + 



2 , *i Xi 



[m~ — mf^ — ml) + w (m~ — m~ — mf) 2 — 4m? m\ 



= 475 GeV. (5.14) 

However searching the endpoint in a tagged tb distribution is not successful. 
This is because there are typically 4 b partons in the final state. When an event has 



3 The events may be double counted in the different regions. 
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Figure 14: The Wbj invariant mass distribution after the W sideband subtraction. The left 
(right) figure is a result of parton (jet) level analysis. 



2 tagged b jets in the final state and one of the jjb system has a mass consistent 
with top mass, the probability that the other b jets is coming from the same gluino 
cascade decay chain is only 1/3. In addition, the probability that both the two b 
tagged jets are satisfy the hemisphere cut pr > 50 GeV and \rj\ < 2.5 is not large. 
One would also loose events by requiring m(bW) ~ m t . 

Instead, we study the Wbj distirbution. The W is selected by the sideband 
method, while one of the other two jets in the same hemisphere is tagged as b. The 
other jet is the highest jet in the hemisphere except for the used jets. When there 
are more than one b tagged jets in a hemisphere, we take the highest px b jet. 

We show the Wbj invariant mass distributions after the W sideband subtraction 
in Fig. O. The left figure shows the result of a parton level analysis. Here the 
parton momenta from cascade decays are taken from HERWIG event record and 
they are processed by the same hemishpere algorithm. We assume 60% tagging 
efficiency for b quark. The total distribution has a clear endpoint of the cascade 
decay, g — > tit — > xfbt, expected at 475 GeV. The contribution of the events where 
at least one gluino decays though the decay chain (a) is shown by a dashed line in 
the figure. A dotted line represents the contribution from the events that contains 
subdominant decay chain (d), which does not show such structure at ~ 475 GeV. 
The other gluino decay modes do not give statistically significant contributions. 

The jet level distribution is shown in Fig. [14] (right). The distribution is smeared 
due to the jet energy resolution. Note that this is a distribution of 4 jet system, 
where two of the jets are from b parton. Nevertheless, the endpoint structure is still 
visible. We fit the fitting function used in Ref. ]7l| to the distribution. The fitted 
endpoint is at 468 ± 7 GeV, which is consistent with the theoretically expected value 
475 GeV. The significant edge proves the existence of the decay g — > tit — > x^tb. 
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5.2 Inclusive Mr and M™ 2 in distributions 

Recently, it is realized that M^ 2 variable is useful to determine the squark and gluino 
masses at the LHC |]32], The M T2 variable is defined for a system of sparticle 
pair production and decay, namely the system with two visible objects p™ and p^f. 
and missing momentum arising from 2 LSPs pV = Pj" SP1 +Pj>2 P2 > as a function of an 
arbitrary test LSP mass m x , as follows: 

M T2 (m x ) = . min [ ma^m? (p™, p?f), m?(p^, #)}] , (5.15) 

_miss i _miss - 1 

PTl +PT2 -PT 

The minimization is taken for the test LSP momenta, p™} ss and p™ 2 Sf \ under the 
constraint p™f s + p™ 2 ss — p'r- The transverse mass, mg , is defined as 

[™?(Pt* , p?f)] 2 = Kf + m 2 + 2(£™£— - p™ ■ p™»), (5.16) 

where mj is is the invariant mass of the "visible object", (m™ 3 ) 2 = (pj is ) 2 and E?i = 
p 2 ^ + m 2 . The kinematical upper bound of the transverse mass is given by 



the mother particle mass if p™ ss = p^. and m x = m^o. Because of this property 
the kinematical upper bound of the Mti variable is given by 

Mr2(w^o) < max{mi, TO2}, (5-17) 

where m\ and m 2 are masses of the initially produced sparticles. Thus, the endpoint 
of the Mt2 distribution relates to the heavy squark mass in our scenario as: 

M^ d (m x o) ~ rriq. (5.18) 

Squark and gluino production events often produce O(10) jets in the final state. 
The central question for the application of Mt2 analysis is how to define the p^f 
in such cases. The "inclusive Mt2 V is defined so that p^ ls is taken as a hemisphere 
momentum in Eq. ( |5.10|) . In addition, a sub-sytem Mt2 is introduced in Refs. |33], 



|34| which is defined as an inclusive M T2 but the highest px jet is removed before 
hemisphere reconstruction. As discussed in the previous section, a quark jet from 
the decay q ^ qg or q ^ qxi tends to be the highest pr jet if 3> nig, which is also 
the case for our scenario (See Eq. ( |3.2p ). The sub-system after removing the highest 
jet tends to be g-g or 'g-Xi system for 'g-q production events. Thus the endpoint of 
the subsystem M T2 should be the gluino mass. 

More systematical approach to observe the gluino mass is using the minimum 
Mt2- The minimum Mt2 {M^ 2 m ) is defined by 

M™ in = . min [M^f{i)] , (5.19) 

where M^ 2 {i) is a generalized sub-system Mt2 defined by removing the i-th high pt 
jet before hemisphere reconstruction. The M™ 2 m has been defined in Ref. |36| for the 
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Figure 15: The inclusive M T 2 distributions at Points A, B, AH1 and AH2. The result of the 
endpoint fit are M^ n 2 d = 1199 ± 23 and 1051 ± 16GeV for Point A and B, respectively. The first 
two generation squark masses are 1150 and 1080 GeV for Points A and B, respectively. 



leading 5 jets of the events pp —>gg — > 4j + 2xi to reduce an effect of the initial state 
radiation (ISR) to a g-g system. Namely, if a gluino pair is produced via a gq — > ggj 
process, this jet may have a large px compared with some of jets from a sparticle 
cascade decay. By using M™ 2 in one can effectively reduces the contamination from 
the ISR jet in g-g production events. 

In this paper, we use M™ 2 in to reduce the additional jet from q — > qg in g-q events 
as well as the effect of the ISR jet in g-g events. Note that for our model points, a 
gluino decay may lead 6 jets in the final state, therefore the chance that those jets 
are soft compared with the ISR is not small. The endpoint of the M™ 2 m is given by 
the gluino mass: 

(M£fK ? )r d ~m ? . (5.20) 

We show the inclusive My 2 distributions of our sample model points in Fig. |15"|. 
Here we choose m x = 100 GeV, and p^ s = p^ mi and p^ s = p^eLi for the Mti 
calculation. 9 In addition to the standard SUSY cut, we require no isolated lepton, 
Py* > 50 GeV and \if el \ < 2.5 for the jet involved in a hemisphere. We also require 
Nino > 1 for Points A and B, and iV*J (600) > 1 for Point AH1 (AH2) to select 

q-g(q) production events. We can optimize this cut from the difference of the p^p 
distributions for b tagged and non-6 tagged jets (See Figs. || and |6p. The cross sections 
of squark productions are small for Points AH1 and AH2 (See Table|^). We use events 
correspond to 20 fb -1 for Points AH1 and AH2, while 5fb _1 for Points A and B. 

The Mt2 distributions at Point A and B have endpoints near the squark mass: 
rriq ~ 1150 and 1080 for Points A and B, respectively. We fit the distributions by a 
simple fitting function 

/(m) = 0(m - M^ d )[ai(m - M^ d ) + b} + 0(M^ d - m)[a 2 {M^ 2 d -m)+b], (5.21) 

9 We require /^(seed) and p 2 (seed) remain in the different hemispheres when we calculate Mj- 2 
and M™ 2 m . The condition seems important to keep the events near the endpoint of M™ 2 m . 
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Figure 16: The distributions of the M T2 (= min{ M$,M$}) at AH1 and AH2. The results 
of the endpoint fit are M™ d = 1530 ± 31 and 1798 ± 21 GeV for Points AH1 and AH2, respectively. 
The first two generation squark masses are rrig ~ 1500 and 1780 GeV for Points AH1 and AH2, 
respectively. 



to see if the endpoints are recovered correctly. We obtain Mf2, d = 1199 ± 23 and 
1051 ± 16 GeV for Points A and B, respectively. They are roughly consistent with 
the input squark masses. 

For Points AH1 and AH2 (especially for AH2), there are a few events in the 
regions where Mr 2 < rrtq. At these model points, a half of the heavy squark mass 
is larger than the gluino mass, m§-/2 > m-g. It is therefore expected that a gluino 
from the squark decay q — > gq is boosted and goes in the opposite direction from the 
direction of the quark jet q. In that case, decay products of the g and the quark jet 
q are not likely to be in the same hemisphere. 

To reconstruct the heavy squark mass, we should separate the final states of 
g-q — > g-g-q events into three parts, two groups of decay products of the gluinos and 
the quark jet q. We adopt the same technique as that used in the sub-system Mt2- 
Namely, we remove the highest pt jet in the event, and the rest of jets and leptons 
are grouped into the each hemisphere. Next we assign the highest pt jet into one 
of the hemispheres and calculate Mr 2 . This gives two Mj> 2 values, and M^ , 
depending on which hemisphere the quark jet is assigned. Finally, we choose the 
smaller M^, M T2 = mm{M$, M^}. 

The distributions of the Mt2 calculated from such procedure at Points AH1 
and AH2 are shown in Fig. Here we adopt the same cuts as in Fig. pl| The 
distributions have robust endpoint structures near the correct squark masses: 1500 
and 1780 GeV for Points AH1 and AH2, respectively. We fit the distributions by the 
fitting function in Eq. flOl)) . The results are M^f = 1530 ± 31 and 1798 ± 21 GeV 
for Points AH1 and AH2, respectively. They are consistent with the input squark 
masses. 
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Figure 17: The Af™ 2 ln distributions at Points A and B. The results of the endpoint fit are 
(MJ? 2 in ) cnd = 715 ± 14 and 524 ± 8GeV for Points A and B, respectively. The gluino masses arc 
nig ~ 697 and 544 GeV for Points A and B, respectively. 



Next, we show the M™ 2 m distributions for Points A and B in Fig. [L7|. Here, 
we require = to reduce the q-q(g) production events. The numbers of events 
corresponds to f £ = 5fb _1 for each model points. We fit the f{x) in Eq. ( 5.21 



to the Mf 2 in distributions. We obtain (Mf 2 in ) end = 715 ± 14 and 524 ± 8 GeV at 
Points A and B, respectively. They are roughly consistent with input gluino masses: 
m-g = 697 and 544 GeV for Points A and B, respectively. 

If one assume GUT relation among the gaugino masses and gaugino dominance 
of the Xi > one can roughly estimate the Xi mass from the gluino mass. Since the tb 
endpoint is a function of m-g, m^± and , one can get the information of from 
the the measured gluino mass and the tb endpoint. 



6. Summary and Conclusion 

In this paper, we investigate the LHC signature of modified universal sfermion mass 
(MUSM) scenario. In the scenario the sfermion mass matrices can be parametrized 
as in Eq. (|2.3| ) at the GUT scale. In this paper, we concentrate the region where 
m ^> m 30 , rrii/2, m H u , H based on the considerations of naturalness and flavor and 
CP constraints. 

In this scenario, gluino decays entirely into the 3rd generation squarks. The 
SUSY events typically have 4 b partons. The fraction of the SUSY events without b 
tagged jets is suppressed even if the b tagging efficiency is 60%. This is a feature of 
models with the mass relation Eq. ( |3.1| ). 

The 1st and 2nd generation squarks are much heavier than gluino and the 3rd 
generation squarks in this scenario. The mass of the heavy squarks mo cannot be 

10 We do not show the distributions for Points AH1 and AH2 because their gluino mass are almost 
the same as that of Point A. 
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arbitrarily large from CCB constraint. The heavy squarks can be observed at the 
LHC if m < 1 — 2 TeV. A quark jet from the decay q — > gq tends to have large px 
(Pt ~ m g/2), and they will be tagged as the excess of the non-6 tagged jets in the 
high p T range relative to the b tagged jets. The excess indicates the mass relation 
Eq.(23). 

The signature similar to that of MUSM scenario may be observed at some 
CMSSM points with \Aq\ > mo ^> mm. For MUSM scenario, gluino decay modes 
are g — ► tit and g — > bib, while gluino decay is dominated by g —>■ tit mode for the 
CMSSM region . This difference may be seen in the \\ ~ > Xi^ channel. If a \\ is 
originated from the bi decay, the number of associated jets is 2 (bb). On the other 
hand if X2 ls originated from the t% decay, it is 6 (bb + 2W(jj)) when both of the W 
bosons decay into jets. We demonstrate that we can discriminate these scenarios by 
investigating the number of jets in the event with 2 leptons with m\i ~ mz- 

In MUSM scenario, the events contain many jets arising from gluino decays. The 
mass reconstructions of the SUSY particles are challenging due to the combinatorial 
background. However, we find that successful reconstructions are possible. The con- 
ventional endpoint analysis with OSSF lepton pair from a cascade decay chain (|5.4| ) 
is useful if rain and \A Q \ are small enough to satisfy conditions (|5.3| ). We demon- 
strate that all sparticle masses arising from the decay chain ( |5.4|) can be measured 
at the LHC in such case. On the other hand, if mi/2 or \Aq\ is large enough, we 
have to use events without leptons for the exclusive analysis. We have succeeded to 
reduce the combinatorial background by searching for jet pair consistent with W in 
the same hemisphere. Especially we can efficiently reconstruct the top quark arising 
from gluino decays into scalar top. Moreover we demonstrate that the endpoint in 
the Wbj distribution can be seen even in the jet level analysis, which indicated the 
tb endpoint of the g — > tit — > xfbt mode. 

The inclusive M^ 2 and M™ 2 m distribution is also useful for the mass determination 
of gluino and the heavy squarks. The squark mass can be measured from the endpoint 
of the M T2 distribution. Moreover we reconstruct the 'g-g and/or g~Xi "sub-system" 
by removing the 2-th high pt jet before hemisphere reconstruction and calculating 
subsystem Mt2 called M T ^(i). The gluino mass can be measured from the endpoint 
of the M™ in (= min^ =1 [Af|i2 b (0] ) distribution. By combining the measured gluino 
mass and tb endpoint, one can get the information of the stop mass only from the 
analysis using jets. Such inclusive analyses may tell us both the mass scales of 
squarks and gluino in the early stage of the LHC such as J £ ~ 5 — 20 fb _1 . 
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A. Masses and Branching Ratios 



Table 4: Left; Mass spectra of sparticles and Higgs bosons for our model points in GeV. Right; 
Branching ratios of sparticles for our model points. Here, u and d denote both the first and the 
second generation up and down type quarks, respectively 
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B. Comparison with SPS benchmark points 



For comparison of MUSM scenario with other scenarios with universal sfermion 
masses, we show the distribution of number of b tagged jet and the distribu- 
tions at several benchmark model points called snowmass points and slopes (SPS) 
pTj|j . Here, the same cut as in Section [3] is adopted. The number of generated events 
correspond to J Cdt = lfb -1 . 



SPS la "Typical point" 
m = 100 GeV, m 1/2 = 250 GeV, A = -100 GeV, tan/? = 10, sgn(/i) = + 
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SPS3 "Coannihilation region" 

m = 90 GeV, m 1/2 = 400 GeV, A Q = GeV, tan (3 = 10, sgn(/i) = + 
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SPS4 "Large tan/3" 

m = 400 GeV, m 1/2 = 300 GeV, A = GeV, tan (3 = 50, sgn(/i) = + 
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SPS5 "Light stop" 
m = 150 GeV, m 1/2 = 300 GeV, A = -1000 GeV, tan = 5, sgn(/i) = + 
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SPS6 "Non-universal gaugino" 
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SPS9 "AMSB scenario" 



m = 400 GeV, m 3/2 = 60 TeV, tan (3 = 10, sgn(/x) = + 



SPS9 




Number of b-tagged jets 



SPS9 

i 

I 



rL 



ran 



200 400 600 800 1000 1200 

(GeV) 

The highest pT (not b— tagged) 



SPS9 



I 200 400 600 800 1000 1200 

(GeV) 

The highest pT (b-tagged) 



- 32 - 



References 

[1] H. P. Nilles, Phys. Rept. 110, 1 (1984). 

[2] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

[3] S. P. Martin, arXiv:hep-ph/9709356. 

[4] I. Hinchliffe, F. E. Paige, M. D. Shapiro, J. Soderqvist and W. Yao, Phys. Rev. D 55 
(1997) 5520 [arXiv:hep-ph/9610544]. 

[5] I. Hinchliffe and F. E. Paige, Phys. Rev. D 61, 095011 (2000) 
[arXiv:hep-ph/9907519]. 

[6] H. Bachacou, I. Hinchliffe and F. E. Paige, Phys. Rev. D 62 (2000) 015009 
[arXiv:hep-ph/9907518]. 

[7] B. C. Allanach, C. G. Lester, M. A. Parker and B. R. Webber, JHEP 0009, 004 
(2000) [arXiv:hep-ph/0007009]. 

[8] J. L. Feng and T. Moroi, Phys. Rev. D 58, 035001 (1998) [arXiv:hep-ph/9712499]. 

[9] I. Hinchliffe and F. E. Paige, Phys. Rev. D 60, 095002 (1999) 
[arXiv:hep-ph/9812233]. 

[10] K. Kawagoe, T. Kobayashi, M. M. Nojiri and A. Ochi, Phys. Rev. D 69, 035003 
(2004) [arXiv:hep-ph/0309031]. 

[11] F. E. Paige and J. D. Wells, arXiv:hep-ph/0001249. 

[12] A. J. Barr, C. G. Lester, M. A. Parker, B. C. Allanach and P. Richardson, JHEP 
0303, 045 (2003) [arXiv:hep-ph/0208214]. 

[13] A. Datta and K. Huitu, Phys. Rev. D 67, 115006 (2003) [arXiv:hep-ph/0211319]. 

[14] S. Asai, T. Moroi, K. Nishihara and T. T. Yanagida, Phys. Lett. B 653, 81 (2007) 
[arXiv:0705.3086 [hep-ph]]. 

[15] S. Asai, T. Moroi and T. T. Yanagida, Phys. Lett. B 664, 185 (2008) 
[arXiv:0802.3725 [hep-ph]]. 

[16] H. Baer, E. K. Park, X. Tata and T. T. Wang, Phys. Lett. B 641, 447 (2006) 
[arXiv:hep-ph/0607085]. 

[17] W. S. Cho, Y. G. Kim, K. Y. Lee, C. B. Park and Y. Shimizu, JHEP 0704, 054 
(2007) [arXiv:hep-ph/0703163]. 

[18] S. Abdullin et al. [CMS Collaboration], J. Phys. G 28 (2002) 469 
[arXiv:hep-ph/9806366]. 

[19] ATLAS Collaboration, "ATLAS detector and physic performance Technical Design 
Report," CERN/LHCC 99-14/15 (1999). 



- 33 - 



[20] B. K. Gjelsten, D. J. 2. Miller and P. Osland, JHEP 0412, 003 (2004) 
[arXiv:hep-ph/0410303]. 

[21] D. J. 2. Miller, P. Osland and A. R. Raklev, JHEP 0603, 034 (2006) 
[arXiv:hep-ph/0510356]. 

[22] B. K. Gjelsten, D. J. 2. Miller, P. Osland and A. R. Raklev, AIP Conf. Proc. 903, 
257 (2007) [arXiv:hep-ph/0611259]. 

[23] M. M. Nojiri, G. Polesello and D. R. Tovey, arXiv:hep-ph/0312317. 

[24] K. Kawagoe, M. M. Nojiri and G. Polesello, Phys. Rev. D 71, 035008 (2005) 
[arXiv:hep-ph/0410160]. 

[25] M. M. Nojiri, G. Polesello and D. R. Tovey, JHEP 0805, 014 (2008) 
[arXiv:0712.2718 [hep-ph]]. 

[26] H. C. Cheng, D. Engelhardt, J. F. Gunion, Z. Han and B. McElrath, Phys. Rev. 
Lett. 100, 252001 (2008) [arXiv:0802.4290 [hep-ph]]. 

[27] C. G. Lester and D. J. Summers, Phys. Lett. B 463, 99 (1999) 
[arXiv:hep-ph/9906349]. 

[28] A. Barr, C. Lester and P. Stephens, J. Phys. G 29, 2343 (2003) 
[arXiv:hep-ph/0304226]. 

[29] W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, Phys. Rev. Lett. 100, 171801 (2008) 
[arXiv:0709.0288 [hep-ph]]. 

[30] A. J. Barr, B. Gripaios and C. G. Lester, JHEP 0802, 014 (2008) [arXiv:071 1.4008 
[hep-ph]]. 

[31] W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, JHEP 0802, 035 (2008) 
[arXiv:071 1.4526 [hep-ph]]. 

[32] M. M. Nojiri, Y. Shimizu, S. Okada and K. Kawagoe, JHEP 0806, 035 (2008) 
[arXiv:0802.2412 [hep-ph]]. 

[33] M. M. Nojiri, K. Sakurai, Y. Shimizu and M. Takeuchi, JHEP 0810, 100 (2008) 
[arXiv:0808.1094 [hep-ph]]. 

[34] M. Burns, K. Kong, K. T. Matchev and M. Park, JHEP 0903, 143 (2009) 
[arXiv:0810.5576 [hep-ph]]. 

[35] W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, arXiv:0810.4853 [hep-ph]. 

[36] J. Alwall, K. Hiramatsu, M. M. Nojiri and Y. Shimizu, arXiv:0905.1201 [hep-ph]. 

[37] S. I. Bityukov and N. V. Krasnikov, Phys. Atom. Nucl. 65, 1341 (2002) [Yad. Fiz. 
65, 1374 (2002)] [arXiv:hep-ph/0102179]. 



- 34 - 



[38] S. Bhattacharya, A. Datta and B. Mukhopadhyaya, Phys. Rev. D 78, 035011 (2008) 
[arXiv:0804.4051 [hep-ph]]; Phys. Rev. D 78, 115018 (2008) [arXiv:0809.2012 
[hep-ph]]. 

[39] R. Barbieri and D. Pappadopulo, arXiv:0906.4546 [hep-ph]. 

[40] J. R. Ellis and D. V. Nanopoulos, Phys. Lett. B 110, 44 (1982); 

[41] R. Barbieri and R. Gatto, Phys. Lett. B 110, 211 (1982). 

[42] J. S. Hagelin, S. Kelley and T. Tanaka, Nucl. Phys. B 415, 293 (1994). 

[43] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477, 321 
(1996) [arXiv:hep-ph/9604387]. 

[44] S. Dimopoulos and G. F. Giudice, Phys. Lett. B 357, 573 (1995) 
[arXiv:hep-ph/9507282]. 

[45] A. Pomarol and D. Tommasini, Nucl. Phys. B 466, 3 (1996) [arXiv:hep-ph/9507462]. 

[46] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 388, 588 (1996) 
[arXiv:hep-ph/9607394] . 

[47] S. G. Kim, N. Maekawa, A. Matsuzaki, K. Sakurai and T. Yoshikawa, Phys. Rev. D 
75, 115008 (2007) [arXiv:hep-ph/0612370]. 

[48] S. G. Kim, N. Maekawa, A. Matsuzaki, K. Sakurai and T. Yoshikawa, Prog. Theor. 
Phys. 121, 49 (2009) [arXiv:0803.4250 [hep-ph]]. 

[49] M. Ishiduki, S. G. Kim, N. Maekawa and K. Sakurai, arXiv:0901.3400 [hep-ph]. 

[50] N. Maekawa, Phys. Lett. B561, 273 (2003)[arXiv:hep-ph/0212141]; Prog. Theor. 
Phys. 112, 639 (2004)[arXiv:hep-ph/0402224]; 

[51] N. Maekawa, T. Yamashita, JHEP 0407, 009 (2004)[arXiv:hep-ph/0404020]. 

[52] J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. B 510, 236 
(2001) [arXiv:hep-ph/0102098]; L. Roszkowski, R. Ruiz de Austri and T. Nihei, 
JHEP 0108, 024 (2001) [arXiv:hep-ph/0106334]; J. R. Ellis, K. A. Olive and 
Y. Santoso, New J. Phys. 4, 32 (2002) [arXiv:hep-ph/0202110]. 

[53] R. Dermisek, H. D. Kim and I. W. Kim, JHEP 0610, 001 (2006) 
[arXiv:hep-ph/0607169]. 

[54] K. Hagiwara, A. D. Martin, D. Nomura and T. Teubner, Phys. Lett. B 649, 173 
(2007) [arXiv:hep-ph/0611102] and references are there in. 

[55] N. Arkani-Hamed and H. Murayama, Phys. Rev. D 56, 6733 (1997) 

[arXiv:hep-ph/9703259]; K. Agashe and M. Graesser, Phys. Rev. D 59, 015007 
(1999) [arXiv:hep-ph/9801446]. 



- 35 - 



[56] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, arXiv:hep-ph/0312045. 

[57] R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH 
Collaboration and and], Phys. Lett. B 565, 61 (2003) [arXiv:hep-ex/0306033]. 

[58] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[59] Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 85, 1 (1991); Phys. 
Lett. B 262, 54 (1991); J. R. Ellis, G. Ridolfi and F. Zwirner, Phys. Lett. B 257, 83 
(1991); Phys. Lett. B 262, 477 (1991); Phys. Lett. B 271, 123 (1991); H. E. Haber 
and R. Hempfhng, Phys. Rev. Lett. 66, 1815 (1991); Phys. Rev. D 48, 4280 (1993) 
[arXiv:hep-ph/9307201]; M. Drees and M. M. Nojiri, Phys. Rev. D 45, 2482 (1992); 
J. R. Espinosa and R. J. Zhang, JHEP 0003, 026 (2000) [arXiv:hep-ph/9912236]; 
M. S. Carena, H. E. Haber, S. Heinemeyer, W. Hollik, C. E. M. Wagner and 
G. Weiglein, Nucl. Phys. B 580, 29 (2000) [arXiv:hep-ph/0001002]. 

[60] G. L. Kane, T. T. Wang, B. D. Nelson and L. T. Wang, Phys. Rev. D 71, 035006 
(2005) [arXiv:hep-ph/0407001]; M. Drees, Phys. Rev. D 71, 115006 (2005) 
[arXiv:hep-ph/0502075]; A. Belyaev, Q. H. Cao, D. Nomura, K. Tobe and 
C. P. Yuan, Phys. Rev. Lett. 100, 061801 (2008) [arXiv:hep-ph/0609079]. 

[61] S. G. Kim, N. Maekawa, A. Matsuzaki, K. Sakurai, A. I. Sanda and T. Yoshikawa, 
Phys. Rev. D 74, 115016 (2006) [arXiv:hep-ph/0609076]; 

[62] M. Asano, S. Matsumoto, M. Senami and H. Sugiyama, Phys. Lett. B 663, 330 
(2008) [arXiv:0711.3950 [hep-ph]]. 

[63] S. G. Kim, N. Maekawa, K. I. Nagao, K. Sakurai and T. Yoshikawa, Phys. Rev. D 
78, 075010 (2008) [arXiv:0804.3084 [hep-ph]]. 

[64] W. S. Hou, Phys. Rev. D 48, 2342 (1993). 

[65] K. S. Babu and C. F. Kolda, Phys. Rev. Lett. 84, 228 (2000) 
[arXiv:hep-ph/9909476]. 

[66] G. Isidori and P. Paradisi, Phys. Lett. B 639, 499 (2006) [arXiv:hep-ph/0605012]; 
J. R. Ellis, S. Heinemeyer, K. A. Olive, A. M. Weber and G. Weiglein, JHEP 0708, 
083 (2007) [arXiv:0706.0652 [hep-ph]]; F. Domingo and U. Ellwanger, JHEP 0712, 
090 (2007) [arXiv:0710.3714 [hep-ph]]. 

[67] G. Corcella et al, JHEP 0101, 010 (2001) [arXiv:hep-ph/0011363]. 

[68] G. Corcella et al, arXiv:hep-ph/0210213. 

[69] E. Richter-Was, arXiv:hep-ph/0207355. 

[70] B. C. Allanach et al, in Proc. of the APS/DPF/DPB Summer Study on the Future 
of Particle Physics (Snowmass 2001) ed. N. Graf, Eur. Phys. J. C 25, 113 (2002) 
[arXiv:hep-ph/0202233]. 



- 36 - 



[71] K. i. Hikasa and M. Kobayashi, Phys. Rev. D 36, 724 (1987). 

[72] M. Carena, A. Freitas and C. E. M. Wagner, JHEP 0810, 109 (2008) 
[arXiv:0808.2298 [hep-ph]]. 

[73] A. Bartl, W. Majerotto and W. Porod, Z. Phys. C 64, 499 (1994) [Erratum-ibid. C 
68, 518 (1995)]. 

[74] U. Chattopadhyay, A. Datta, A. Datta, A. Datta and D. P. Roy, Phys. Lett. B 493, 
127 (2000) [arXiv:hep-ph/0008228]. 

[75] J. Hisano, K. Kawagoe, R. Kitano and M. M. Nojiri, Phys. Rev. D 66, 115004 

(2002) [arXiv:hep-ph/0204078]; J. Hisano, K. Kawagoe and M. M. Nojiri, Phys. Rev. 
D 68, 035007 (2003) [arXiv:hep-ph/0304214]. 

[76] F. Moortgat and L. Pape, CMS Physics TDR, Vol. II, Report No. 
CERN-LHCC-2006, Chap. 13.4, p410 

[77] S. Matsumoto, M. M. Nojiri and D. Nomura, Phys. Rev. D 75 (2007) 055006 
[arXiv:hep-ph/0612249]. 

[78] J. Hubisz, J. Lykken, M. Pierini and M. Spiropulu, arXiv: 0805. 2398 [hep-ph]. 



-37- 



